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ABSTRACT 

We compute the escape of ionizing radiation from galaxies in the redshift interval ^ = 4 — 10, i.e., 
during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully 
cosmological simulations. The simulations invoke early, energetic feedback, and the galaxies evolve 
into a realistic population at z = 0. Our galaxies cover nearly four orders of magnitude in masses 
(lO'''® — lO^"'^'^ M0) and more than five orders in star formation rates (10~^'^ — 10^'^ yr~^), and we 
include an approximate treatment of dust absorption. We show that the source-averaged Lyman-limit 
escape fraction at z = 10.4 is close to 80% declining monotonically with time as more massive objects 
build up at lower redshifts. Although the amount of dust absorption is uncertain to 1 — 1.5 dex, it is 
tightly correlated with metallicity; we find that dust is unlikely to significantly impact the observed 
UV output. These results support reionization by stellar radiation from low-luminosity dwarf galaxies 
and are also compatible with Lyman continuum observations and theoretical predictions at z 3 — 4. 

Subject headings: galaxies: dwarf — galaxies: formation — methods: numerical — stars: formation 



1. INTRODUCTION 

The de clining number density of quasars at redshift 
z ~> 2.5 (|Fan et al.ll200lD points to stellar radiation as 
a likely source of cosmic hydrogen reionization. One 
of the key unknown parameters of reionization mod- 
els is the escape fraction Jcsc of ionizing photons from 
high-redshift galaxies. Observationally, f^sc is much 
easier to estimate at lower redshifts. In local star- 
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factors with hydrogen column densities of ~> 10^^cm~^ 
and escape of ionizing radiation through relatively few 
transparent holes in the interstellar medium (ISM) . Star- 
bursts in the disk of our own Galaxy feature /esc ^ 6% 
normal to the disk (jBland-Hawthorn fc MalonevI 119991 
[2001). a value which might be a factor of 2 — 3 lower 
when averaged over all directions. 

There is strong evidence that the Lyman continuum 
(LyC) escape fraction evolves with redshift. A quan- 
tity often measured in high-redshift observations is the 
relative escape fraction fl°]. — /osc//csc,i500 defined as 
the ratio of the escape fraction of Lyman-limit pho- 
tons /esc to that of I5OOA photons. The conversion be- 
tween /esc and the relative escape fractions depends on 
the amount of dust attenuati on at 1500A with in each 
galaxy and is often uncertain (jinoue et al.l 120051 ). Large 
observed samples of z ~ 1 starburst galax ies point to 
fairly low /^"^ < 8 - 10% (ISiana et al.l(2007l) : however, 
the ionizing output changes toward higher redshift. At 
2^3, from direct detection of LyC emission from 2 of 
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the 14 studied star-forming galaxies fShaplev et al.l (|2OO60 
put constraints on the relative escape fraction around 
fill ~ 14%. More recently, Iwat a et al.l (j2009) expanded 
the search to 198 galaxies at z w 3.1, detecting LyC ra- 
diation from seven Lyman break galaxies (LBGs) and 
10 Lyman a emitters (LAEs), estimating that prior to 
accounting for intergalactic medium (IGM) attenuation 
the average relative escape fraction is fUl = 61%. Tak- 
ing into account IGM attenuation with a median LyC 
opacity tigm = 0.59 results in fl°l — 110%, which in 
their estimate translates into the average absolute es- 
cape fraction f^sc = 26%, and a lower limit of /esc = 15% 
prior to accounting for IGM attenuation. Furthermore, 
from a comparison of direct observations of the LyC 
from galaxies and indirect estimates f r om in tergalactic 
medium (IGM) ionization, llnoue "eFall (|2OO60 found ev- 
idence for redshift evolution of the average /esc from 
1 - 2% at z = 2 to - 10% at z ^ 3.6. 

Some of the first theoretical estimates of /esc were 
derived from studies of gas distribution in the Milky 
Way and from anal ytical models of galaxy formation. 
IDove h Shuil (fl99l found /esc ~ 7% for OB associa- 
tions in the Galactic disk. Additional absorption comes 
from the shells of expanding superbubbles which further 
become Rayleigh- Taylor unstable, making the problem of 
calculating /esc dependent o n time and on details of star 
formation (I Dove et al.ll2000l ) . Analytical estimates of /esc 
in the cosmological context allowed us to study its de- 
pendence on multiple parameters, such as redshift, dark 
matter halo mass, baryonic fractio n, star formation effi- 
ciency, and stellar ionizing output (jRicotti fc Shulll[2000l : 
Wood & Loeb 2000). 

More accurate theoretical constraints on /esc can be 
derived from cosmological hydrodynamical models in- 
corporating star formation and transfer of ionizing ra- 
diation in the surrounding gas. At very high redshifts 
(z ~ 20) estimates show that the escape fractions refiect 
the large masses of Population III s tars and can easily 
approach unity (Whalen et a l. 2004; Alvarez et al.l[2006.' ) . 
In the observable, post-reionization universe, there have 
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been several attempts to compute /esc from cosmologi- 
cal galaxy formation models. In our earl ier simulations 
azoumov fc Sommer-Larsenl2006l[200l we found that 
the source-averaged /esc monotonically declines with red- 
shift, from - 6 - 10% at ^ = 3.6 to ~ 1 - 2% at z = 2.39 
in lin e with the observational findings of llnoue et al.l 
(120061). 

iGnedin et al.l ()2008[ ) used a high-resolution simula- 
tion of a Lagrangian region of a 6/i~^ Mpc cosmolog- 
ical model. The selected region corresponds to five 
virial radii of a Milky Way-sized galaxy and includes the 
Milky Way progenitor and several dozen smaller galax- 
ies spanning two orders of magnitude in masses and star 
formation (SF) rates. Analyzing these galaxies at in - 
termediate (z = 3 — 5) redshifts, IGnedin et al.l (|2008 h 
found that the average escape fractions do not depend 
on redshift and tend to correlate positively with the 
SF rate dropping below 10~'^ for galaxies with SF rates 
less than ~ 1 Mq yr~^ and dark matter masses below 
~ 10^°-^ Mq. Such low /esc are explained by the fact that 
these galaxies have a relatively small gravitational pull 
and achieve pressure and density necessary for SF only 
near the midplane of a thick neutral disk. The resulting 
SF regions therefore tend to be embedded d e ep ins ide the 
opaque disk. Consequently, IGnedin et al.l (|2008| ) point 
out that high-redshift galaxies might be quite inefficient 
in emitting ionizing radiation into the IGM. 

It is widely agreed that in the absence of energetic 
supernovae-driven outflows the increased disk density 
at higher redshifts z ~ 1 will result in ver y low 
/esc ^ 1% for smooth disks (jWood fc Loebll200flh . On 
the other hand, simulated galactic disks in cold dark 
matter (CDM) cosmologies suffer the excessive loss of an- 
gular momentum resulting in galaxies th at are too s mall 
compared to observations f|Fall fc Efstathioul Il980[ ). A 
plausible solution to this angular momentum problem 
is heating via feedback from stars and supernovae that 
would impede rapid gas cooling and collapse of small 
protogalactic gas clouds at high redshifts, allowing gas 
to preserve a larger fracti on of its angular momentum 
as it settles into the disk. iSommer-Larsen et al.l (|2003D 
showed that a population of realistic disk, lenticular, and 
elliptical galaxies with almost the correct amount of an- 
gular momentum can be obtained self-consistently in cos- 
mological simulations by invoking energetic stellar feed- 
back at early times, either through a much higher SF 
efficiency than in present-day galaxies, a more top-heavy 
initial mass function (IMF), or both. The affected dwarf 
galaxies would not have time to settl e into a disk in which 
feedb ack would be rather inefficient (|Mac Low fc Ferraral 
Il999| ). Instead, they maintain a more spherical or irregu- 
lar shape resulting in a much stronger impact from stellar 
feedback. 

Recently, using very high resolution (0.1 pc) coupled 
radiation hydrodynamics simulations of d warf galaxies 
in the mass range 10^-^ - 10^-^ M© at z = 8. 1 Wise fc CenI 
([2009) showed that the time-averaged /esc can reach up 
to 80% in halos above 10* M© with a top-heavy IMF, i.e., 
in active star-forming galaxies dominated by the atomic 
line cooling. I t is not clear to wh at extent /esc depends 
on the IMF as I Wise fc CctI (|2009f ) kept the SN feedback 
strength equal in their normal and top-heavy IMF calcu- 
lations varying only the ionizing photon production rate. 



On the other hand, lower mass (< 10^'^ M©), molecular 
cooling galaxies in their simulations have much lower SF 
efficiency and do not contribute a significant flux to reion- 
ize the universe. I rrespective of the mass of the galaxy, 
iWise fc Ce n (200^ found that /esc can vary by up to an 
order of magnitude on the timescale of a few Myrs which 
is the dynamica l time of a star-forming molecular cloud. 

lYaiima et al.l ()2009f) also found very large escape 
fractions of order 20 — 60% at z = 3.7 — 7, an- 
alyzing a very high resolution (1024^^) hydrodynami- 
cal model of a supernov a-dominated primordial galaxy 
dM ori fc Umemural l2006') . However, their model is a sin- 
gle isolated galaxy taken out of the cosmological context. 

In thi s paper, we t a ke an approach complimentary to 
that of IWise fc CenI (|2009[ ). Even though we do not 
have their ultra-high resolution and we solve the radia- 
tive transfer and ionization equations separately from 
the equations of hydrodynamics, we examine the es- 
cape of ionizing photons from a set of high-resolution 
reionization-epoch galaxies that evolve into a realistic 
population of galaxies at z = 0. Specifically, a two- 
mode SF prescription is adopted, with both the SF ef- 
ficiency of the "early" mode and the feedback strength 
calib rated to (1) overcome the an gular momentum prob- 
lem (H ommer-Larsen et al.l |2003[) . and to (2) produce 
a reason able cosmic enrichment history that fits obser- 
vations (jSommer-Larsen fc Fvnbol 120081) . This efficient 
high-z feedback drastically changes the appearance of 
star-forming galaxies during the reionization era leading 
to much higher /esc than, e.g., in z = 3 — 5 dwarf galaxies 
(IGnedin et al.ll2008y 

This paper is organized as follows. First, we describe 
our galaxy formation models and the technique to com- 
pute /esc in 5J2J We present the results of our study in- 
cluding the effect of dust absorption in ij3l and summa- 
rize our findings in []4l 

2. MODEL 

2.1. Simulated galaxies 

We use an impr oved version of the TreeSPH code 
(jSommer-Larsen et a l. 2003) to rerun nine galaxies se- 
lected from a larger cosmological simulation. Table[l]lists 
the following model parameters for each galaxy: physi- 
cal box size L, circular velocity Vc sA, z — 0, mass of 
an SPH and a dark matter particle, and gravitational 
smoothing lengths for SPH and dark matter particles. 
The code employs the "conserva tive entropy" formula- 
tion (jSpringel fc Hernquistl 120021 ). includes star forma- 
tion by conversion of SPH particles into star particles, 
and accounts for non-instantaneous recycling of gas and 
heavy elements via type II and type la SNe and via stellar 
mass loss. Chemical evolution is traced using 10 sepa- 
rate elements (H, He, C, N, O, Mg, Si, S, Ca, and Fe). 
All galaxies in Table [1] were computed with the Salpeter 
IMF which produces a more reasonable cosmic enrich- 
ment history than the lower m etal yield Kroupa IMF 
(jSommer-Larsen fc Fvnboll2008l ). Each galaxy was mod- 
eled with at least 10^ particles. By z = these galaxies 
evolve into systems with a wide range of masses, from a 
dwarf galaxy to a super Milky Way sized system, char- 
acteristic of the "field" . In terms of the luminosity func- 
tion, they span a range from super-L* to about 6 — 7 
magnitudes below L* . For analysis we use outputs at 
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TABLE 1 
List of Modeled Galaxies. 



Galaxy 


Box Size 


Vc(z = 0) 


"Isph 




^sph 






kpc 


km s^^ 


10^ Mq 


IO^Mq 


kpc 


kpc 


S29 


250 


205 


1.4 


8.1 


0.29 


0.52 


S33 


250 


180 


1.4 


8.1 


0.29 


0.52 


S33sc 


250 


300 


8.3 


47 


0.53 


0.94 


S41 


250 


150 


1.4 


8.1 


0.29 


0.52 


S87 


100 


132 


0.18 


1.0 


0.15 


0.26 


S108 


100 


131 


0.18 


1.0 


0.15 


0.26 


SlOSsc 


80 


35 


0.022 


0.13 


0.07 


0.13 


S115 


100 


125 


0.18 


1.0 


0.15 


0.26 


SllSsc 


80 


50 


0.038 


0.22 


0.09 


0.16 



z — 10.4, 8.2, 6.7, 5.7, 5.0, and 4.4 which are separated 
by 200 Myrs. Fig. [T] shows the evolution of density dis- 
tribution in five galaxies that cover a range of masses, 
and Fig. [2] displays in addition the distributions of stel- 
lar particles and heavy elements in galaxy S33. 

2.2. Method 

We model SF by creating discrete star "particles" , each 
representing a population of stars born at almost the 
same time in accordance with the Salpeter IMF. The 
stellar UV luminosity is determined using the p opula- 
tion synthesis package Starburst 1999 (Leithere r et al.l 
Il999( ) with continuous SF distributed among all stars 
younger than 34 Myrs. We compute the effect of this 
radiation on the surrounding gas with radiative trans- 
fer on top of a nested set of grids containing three- 
dimensional distributions of physical variables. To cre- 
ate such a set from the SPH simulation datasets, we 
projected each galaxy onto a 128'^ uniform grid, di- 
vided every base grid cell which contains more than 
-^max — 10 SPH particles, and continued this process 
of subdivision recursively so that no cell contains more 
than iVmax gas particles. For the actual transport of 
ionizing radiation around star-forming regions, we em- 
ploy the ray-tracing scheme with adaptively splitting ra- 
dial rays around point sourc es JAbel fc Wandeltl [20021: 
iRazoumov &: Sommer-Larsenll2006[ ). At each source, we 
start with 12 isot ropic rays using HEALPix discretiza- 
tion on the sphere (jGorski et al. 2002) which further split 
into leaf rays as we move away from the source or enter a 
spatially refined region. The local ray angular refinement 
level riray is increased in integer increments to satisfy the 
condition 

^max(^ray) ^ ^5 (-^) 

where r.^ is the distance between the source i and the cell 
j of the grid refinement level Ij in which the photoreac- 
tion number and energy rates are computed, and the r^ax 
is a monotonically increasing function chosen to ensure 
that at least three ray segments connect every cell in the 
volume to the each source. The photoreaction rates are 
then used iteratively to compute ionization equilibrium. 

In this paper we will work only with the absolute es- 
cape fractions defined as a fraction of photons emitted 
by the source at a given energy that reaches the virial 
radius of the host galaxy. At large redshifts, star forma- 
tion is not confined to a few massive galaxies; instead, it 
takes place in a large number of smaller galaxies scattered 
throughout the volume, and in many cases the concept of 



a host halo becomes quite ambiguous. To overcome this 
problem, for each star particle we calculate its individual 
radius of a sphere containing 200 times the average cos- 
mic critical density at that redshift which is effectively 
a local virial radius of that particle. In calculating the 
ionization balance, we trace all rays to the edges of the 
computational volume; however, fesc for each source is 
evaluated at its individual rvir- To isolate the contribu- 
tion to the UVB of sources outside largest halos, we also 
identify the most massive halo in each volume and cal- 
culate source-averaged /esc at its ryir including only star 
particles inside this radius. 

3. RESULTS 

Fig.|3]shows source-averaged Lyman-limit escape frac- 
tions of all galaxies without invoking the effects of dust 
which will be covered in < j3.11 computed both at the 
virial radius of the most massive galaxy and at indi- 
vidual Tvir of each source. There is a clear trend of a 
larger fesc towards higher redshifts. Note that at z < 4.4 
the majority of our simulated galaxies reproduce the ob- 
served gradual decline to fesc 1 — 2% at z = 2.4 
(jRazoumov fc Sommer-Larsenll2007l ) . Exceptions are the 
more massive {vc{z = 0) > 200kms~^) systems S29, 
S33 and S33sc which form an early large disk and reach 
/esc < 10% before z = 4. A very interesting result is the 
large escape fractions of near unity toward s higher red- 
shifts z ^ 8—10 matching the findings of iWise fc CenI 
12009) in globally turbulent dwarf 10** - 10^° Mg galax- 
ies at z = 8. The escape fractions show an approximate 
inverse correlation with both the halo mass and the SF 
rates (Fig. [Hand Fig. [5]), coinciding with the gradual in- 
crease in galaxy mass and the decline in SF efficiency 
with time. Note that these systems are all actively star- 
forming galaxies; however, their properties evolve con- 
siderably from z = 10.4 to 4.4, as they move from the 
early efficient, self-propagating SF in turbulent, metal- 
poor dwarf galaxies to the late inefficient mode of SF 
from pre-enriched gas s upplied to heavy galactic di sks 
by smooth cooling fiows ijSommer-Larsen et al.|[2003l ). 

Our lower mass 10^ — 10^ Mq) galaxies are still 
massive enough so that even at high SF rates feedback is 
not efficient a t expelling most of the gas from the galaxy 
(jWise fc Cen| [2b09) but plays an important role in shap- 
ing small-scale fluctuations in the ISM, leading to high 
/esc- This can be seen in all our low- mass, high-SF rate 
galaxies all of which have fairly large /osc- On the other 
hand, there is a significant spread in properties of indi- 
vidual galaxies as the escape fraction cannot be derived 
exclusively from the halo mass and the SF rate as demon- 
strated in the lower panel in Fig. [5] 

LyC radiation from stars in secondary/satellite galax- 
ies is an important contributor to the cosmic ionizing 
background during the early stages of galaxy assembly 
as evidenced by the difference between the two types of 
escape fractions in Fig.[3|(open and filled circles). On the 
other hand, in massive galaxies the overwhelming major- 
ity of star- forming regions are found inside the disks, and 
the two escape fractions become almost identical. 

In Fig. [SI we use individual star-forming regions to 
derive a redshift-dependent fit /esc(Afhaio)- Our fitting 
function is based on a simple toy model of the optical 
depth T to the center of a uniform neutral gas sphere 
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Fig. 1.— Projected gas density in galaxies S29, S33, S33sc, S108, SlOSsc (top to bottom), at ^ = 10.4, 8.2, 6.7, 5.7, 5.0, and 4.4 (left 
to right). The side of the image in each row is 125, 125, 250, 80, 80 physical kpc (top to bottom). The color scale is logarithmic from 
10~^ gcm^ to 10^^ gcm'^ in all images. 



TABLE 2 

/osc(Mhalo) FITS. 



as 



10.4 
8.2 
6.7 
5.7 
5.0 
4.4 



2.78 X 10'^ 
1.30 X 10-2 
5.18 X 10-3 
3.42 X 10-3 
6.68 X 10-5 
4.44 X 10-5 



0.105 
0.179 
0.244 
0.262 
0.431 
0.454 



inside the virial radius r oc (1 + z)'^M^[f^. We account 
for deviations from a uniform sphere - clumpy gas, non- 
uniform radial distribution, partial ionization - by intro- 
ducing a second parameter with /osc = 6xp(— aM^g^j^), 
where a and [3 are functions of redshift (Tabled]). 

3.1. Dust absorption 

We follow the prescription of iGnedin et al.l (|2008[ ) to 
include an approximate treatment of dust absorption 
based on the dust extinction curve for the Small Mag- 
ellanic Cloud (SMC). We use oxygen abundance as a 
proxy for metallicity, and assume the SMC metallicity 
to be 20% solar. The dust optical depth can be written 



Tdust(A) — ?path'T-HO'o,SMc/(A) 



0.2 



(2) 



where Zpath is the physical photon path length, nn is 
the local neutral or total hydrogen number density, for 
models with complete dust sublimation and no sublima- 
tion, respectively, (To,smc = 10~^^cm^, Z/Zq is oxygen 
abundance relative to solar, ai id /(A) is the parametric 
wavelength-dependent fit from' Gnedin et all ()2008l ). 

Fig. [7] shows our results for the two extreme cases, 
with and without dust sublimation, for two large and one 
dwarf galaxies. In the sublimation model stellar photons 
destroy dust, and therefore with sufficiently large escape 
fractions the effect of dust on ionizing radiation is nil. 
Only when the escape fraction is small due to accumula- 
tion of neutral hydrogen along the line of sight, we can 
see some small effect of dust in the sublimation model, 
e.g. in the massive galaxy S33sc at z = 5, where the 
relative effect of dust is below 0.2%. Its impact is more 
noticeable in the no-sublimation model. The true escape 
fractions probably lie in between these values, as some 
amount of dust is likely to be destroyed by ionizing radi- 
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Fig. 2. — Projected gas density (top panel), stellar particle distribution (center panel), and projected mass-weighted oxygen abundance 
[O/H] (solar units) in galaxy S33 at z = 10.4, 5.7, and 4.4 (left to right). Each image is 125 physical kpc on a side. Both color scales are 
to the base-10 logarithm, and column density in the top panel is measured in gcm~^. 

tion is proportional to the total gas dens i ty, sim ilar to 
the no-sublimation model of iGnedin et al.l ()2008[ ). but is 
approximately 25 time s weaker than the latter. Despite 
such low cross-sections. lYajima et all (2009) find that the 
escape fractions are significantly regulated by interstellar 
dust in their isolated supernova-dominated galaxy. The 
reason for such a large effect is their relatively strong 
metal enrichment producing on average solar metallicity 
ISM at z = 3.7 (Fig. [8]) which in addition is very irregular 
and moderately optically thin for ionizing photons. 

Observations of UV-selected star-forming galaxies at 
z ~ 2 point to metallic ities in the range 0.1-0.5 dex below 
solar ()Erb et al]|2006[ ). If star formation peaks at z ~ 
2 — 3, metallicities should be substantially lower in the 
interval z ^ 10.4 — 4.4 we study here. All models in our 
set feature a gradual metal enrichment such that most 
of the massive galaxies have [O/H] ~ 0.1 — 0.5 at z = 2 
(Fig. [8]). Since the no-sublimation models put an upper 
limit on dust absorption (Fig. [T]), we conclude that its 
effect is unlikely to significantly impact the output of 
ionizing photons in z J> 4.4 galaxies, except for very rare, 



ation. Similar to IGnedin et al.l ()2008[ ). we find that the 
overall effect of dust is not very large, reducing the escape 
fractions by few percent at most, explained by the sim- 
ple fact that photoionization of neutral hydrogen domi- 
nates the optical depths. Large escape fractions are only 
observed at low hydrogen column densities, where dust 
absorption is usually small. Its effect on the escape of 
ionizing radiation is more noticeable in massive galaxies 
at lower redshifts in which the buildup of heavy elements 
leads to relatively high dust concentrations (Fig. [7]). 

It is worth pointing out that the magnitude of dust ab- 
sorption is highly uncertain a nd can range significantly 
in different studies. Recently, lYaiima "etall (2009) sug- 
gested another simple model of calculating Tdust from 
chemical composition under the following assumptions: 
(1) absorption is caused by dust grains ranging in ra- 
dius from 0.1 /im to 10 ^m; (2) the size distribution of 
dust grains scales as r~^^'^; (3) the density inside dust 
grains is constant at 3gcm~'^; (4) the diffraction effects 
are negligible since ra > 912A ; and (5) the mass den- 
sity of dust is O.OlZ/Z0Pgas- The resulting dust absorp- 
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Fig. 3. — Source-averaged Lyman-limit escape fractions com- 
puted at rvir of the most massive halo in each volume (open circles) 
and at individual r^ir of each star particle (filled circles) vs. red- 
shift for the nine galaxies. The effects of dust are not included in 
this plot. 
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Fig. 4. — Dark matter halo masses Mj^aio (filled squares) and star 
formation rates (star symbols) vs. redshift for the nine galaxies. 

unusually massive systems. 

3.2. Resolution effects 

For the numerical convergence study, we use an ad- 
ditional model K33 to compare the escape fractions at 
z = 5.85 obtained at standard (K33-64) and 8 times the 
mass resolution (K33-512), while keeping the strength of 
stellar feedback and the ionizing luminosity per unit stel- 
lar mass constant. At these two resolutions, this galaxy 
contains 755 star particles of mass 1.42 x 10^ M0, and 
8447 star particles of mass 1.78 x 10'' M0, respectively. 

We find the Lyman-limit f^sc = (0.142,0.118,0.101) 
for K33-64, and (0.215, 0.173, 0.143) for K33-512, at 
r = (0.5,1,2) X Tvir, respectively (Fig. [9]). The star 
formation rate in the higher resolution model is slightly 
larger, 4.4 yr~^ versus 3.2 M© yr~^, which could ex- 
plain the observed differenc e in /esc, since the two quan- 
tities should be coupled (Wise fc Cenll2009l) . In addition. 
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Fig. 5. — Angular averaged escape fraction as a function of the 
host halo mass (top) and the SF rate (center). Bottom: escape 
fractions in the log(SFrate) — log(Miiaio) plane, with the radius 
of each circle proportional to /^g^. The insert box in the bottom 
panel shows /esc = 10~^, 10~^, and 1. Color in all three panels 
represents redshift with the rainbow palette, from 2 = 10.4 (violet) 
to z = 4.4 (red). 

we see the increased porosity of the ISM in the higher res- 
olution model, with more transparent channels between 
dense clumps through which radiation can escape. To be 
conservative, the escape fractions presented in this work 
should probably be considered as lower limits. 

4. DISCUSSION AND CONCLUSIONS 

We have coupled high-resolution galaxy formation 
models with point-source radiative transfer to compute 
the escape fractions of LyC photons from galaxies at 
z = 10.4, 8.2, 6.7, 5.7, 5.0, and 4.4. We confirm very 
lar ge escape fra c tions of near unity at z 8 — 10 found 
bv iWise fc CenI (l2009l) in dwarf 10« - lO^o M© galaxies, 
conducive to efficient stellar reionization. The difference 
between this result and much lower escape fractions in 
-Mhaio ^ 8 X 10^ M© galaxies with similar SF rates at 
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Fig. 6. — The Lyman-limit escape fractions of all star-forming regions (black dots) in all modeled galaxies vs. the host halo mass A/(rvir, z) 
at six redshifts, computed at individual rvir of each star particle. The blue circles and the error bars show /esc averaged inside each mass 
bin, along with standard deviations. The lines show redshift-dependent two- parameter exponential fits (see the text for details). 



z = 3 - 4 (jGnedin et al.ll2008[ ) reflects the model differ- 
ences between reionization-epoch star-forming galaxies 
and lower-redshift dwarf galaxies with a relatively low SF 
efficiency. There are several physical factors that could 
account for very different SF efficiencies in z ~ 6 — 10 
galaxies. These systems are likely to have very low metal- 
licities fo r which the stre ngth of stellar winds is greatly 
reduced (jKudritzkil l2002t ). producing less local disrup- 
tion in star-forming clouds, and leading to higher SF ef- 
ficiency. Lower mctallicities correspond to less efficient 
cooling; however, this effect can be offset by the fact that 
- at a fixed halo mass ~ higher redshift galaxies are more 
compact resulting in more efficient cooling. In addition, 
there are theoretical expec tations of a more top-heavy 
IMF at high redshi fts (A bel et all l2002t iBromm et all 
[200l iPadoan fc Nordlundl l2002t ) . A larger fraction of 
massive stars is likely to yield more efficient feedback on 
galactic scales, producing a larger number of transparent 
channels in the ISM through which ionizing radiation can 
escape the galaxy. 

All our systems are proto-spiral galaxies residing in a 
"field" environment. We did not consider more crowded 
proto-elliptical environments in which evolution is accel- 
erated due to the higher overdensity. In such galaxies the 
UV escape fractions could drop quicker with time than 
what is seen in our current study. In addition, higher 
metallicity might amplify the effect of dust although rel- 
atively few UV photons would escape from such envi- 
ronments. We will address more crowded systems in a 
future study. 

We did not compute the effect of ionizing radiation on 
the hydrodynamical flow which is instead shaped by the 



thermal feedback energy. Would our results change if 
we used a coupled radiation-hydrodynamics approach? 
Ionizing radiation heats up the gas which can then ex- 
pand and leave the star-forming region. iGnedin et al.l 
12008) found overall that coupling of radiative transfer 
and hydrodynamics does not produce a lar ge change in 
.ff,sr- For the dwarf galaxies considered by IWise fc CenI 
(I2009D . they reported fast variations of /esc by up to an 
order of magnitude on the timescale of a few Myrs which 
is the dynamical timescale of a star-forming molecular 
cloud. These changes in /esc are attributed to gas ex- 
pulsion by winds and supernovae creating transparent 
channels for ionizing photons. Since feedback by winds 
and supernovae is already a key part of our hydrody- 
namical calculations, we argue that our postprocessing 
of fixed time snapshots of galaxies with radiative trans- 
fer with the assumption of ionization equilibrium gives 
a fairly accurate estimate of the escape fractions. On 
the other hand, we cannot compute /esc variations on 
short timescales of few Myrs, and therefore we cannot 
tell exactly how representative our computed /esc of the 
time-average escape fractions, something we are hoping 
to address in the future. 

The other possible source of uncertainty in computing 
/esc is gas clumping on sub-resolution scales. Since for 
fixed cloud mass and ionization the optical depth scales 
as with the cloud radius i?, and recombinations are 
more efficient at higher densities, the true /esc could be 
smaller than our current estimates. On the other hand, 
star-forming clouds are shaped by feedback, and the sim- 
ple T (X scaling might not be valid. Moreover, higher 
grid resolution would yield a more porous ISM, with a 
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Fig. 9. — Distribution of star particles by their Lyman-limit /esc 
at individual particle's rvir at two different numerical resolutions 
at z = 5.85. The dashed lines show source-averaged /esc- 



Fig. 7. — Ratio of Lyman-limit escape fractions with complete 
dust sublimation (open diamonds) and no sublimation (filled dia- 
monds) to the escape fractions without dust for two large galaxies 
(top two panels) and a dwarf galaxy (bottom panel), as a function 
of redshift. The solid lines correspond to the ratio of unity. 
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Fig. 8. — Solid, long-dashed and dash-dotted lines show mass- 
weighted oxygen abundance (relative to solar) in galaxies S33, 
S33sc, and S108sc as functions of redshift. The light-shaded area 
represents the range of [O/H] containing 90% of the gas mass in 
S33, i.e., 5% of the gas has metal enrichment above and below 
this region. The dark-shaded area sh ows the range of fO/ H] in 
the supernova-dominated galaxy from lYaiima et al.l II2009I ) , and 
the dotted horizont al line corresponds to metallicity adopted in 
IGnedin et all pOOSi ). 



strong network of transparent channels through which 
radiation can escape. 
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